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ABSTRACT 

We have generated an extensive new suite of massive giant planet atmosphere models and used it to 
obtain fits to photometric data for the planets HR 8799b, c, and d. We consider a wide range of cloudy 
and cloud-free models. The cloudy models incorporate different geometrical and optical thicknesses, 
modal particle sizes, and metallicitics. For each planet and set of cloud parameters, we explore grids 
in gravity and effective temperature, with which we determine constraints on the planet's mass and 
age. Our new models yield statistically significant fits to the data, and conclusively confirm that the 
HR 8799 planets have much thicker clouds than those required to explain data for typical L and T 
dwarfs. Both models with 1) physically thick forsterite clouds and a QO-fim modal particle size and 2) 
clouds made of 1 ^tm-sized pure iron droplets and 1% supersaturation fit the data. Current data are 
insufficient to accurately constrain the microscopic cloud properties, such as composition and particle 
size. The range of best-estimated masses for HR 8799b, HR 8799c, and HR 8799d conservatively span 
2-12 Mj, 6-13 Mj, and 3-11 Mj, respectively and imply coeval ages between ~10 and ~150 Myr, 
consistent with previously reported stellar age. The best-fit temperatures and gravities are slightly 
lower than values obtained by Currie et al. (2011) using even thicker cloud models. Finally, we use 
these models to predict the near-to-mid IR colors of soon-to-be imaged planets. Our models predict 
that planet-mass objects follow a locus in some near-to-mid IR color-magnitude diagrams that is 
clearly separable from the standard L/T dwarf locus for field brown dwarfs. 

Subject headings: planetary systems — planets and satellites: general — planets and satellites: indi- 
vidual (HR 8799b, HR 8799c, HR 8799d) 



1. INTRODUCTION 

Advances in high-contrast imaging have made it pos- 
sible to image young giant planets at wide orbital sep- 
arations around host stars. Major efforts are currently 
underway targeting progressively higher contrast levels, 
e.g., GPI (Macintosh et al. 2006), SPHERE (Beuzit et al. 
2008), Project 1640 (Hinkley et al. 2008,2011), SEEDS 
(Tamura 2009), and NICI (Liu et al. 2010). Current 
instruments are capable of detecting contrast ratios in 
the near infrared down to ~ 10~ 5 (Marois et al. 2010), 
leading to a number of planet detections in recent years. 
The initial discovery of a giant planet around the low- 
mass cool dwarf 2M1207 (Chauvin et al. 2004 ) was fol- 
lowed by more recent discoveries of giant planets around 
massive stellar hosts HR 8799 (Marois et al. 2008,2010), 
Formalhaut (Kalas et al. 2008), and j3 Pic (Lagrange 
et al. 2009). Such detections provide the exciting op- 
portunity of obtaining planetary spectra (Marois et al. 
2008; Lafreniere et al. 2009; Hinz et al. 2010; Janson 
et al. 2010; Bowler et al. 2010) and characterizing their 
atmospheres, in particular their temperature structures, 
compositions, and gravities. Furthermore, the planetary 
radii, masses, and ages can be derived from the atmo- 
spheric parameters (Burrows et al. 2001; Burrows, 2005). 

Extrasolar planets detected with high-contrast imaging 
occupy a unique region in orbital phase space. Planetary 
spectra can also be observed for transiting planets. How- 
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ever, transit surveys for planet detection are sensitive to 
small orbital separations much less than an astronomical 
unit (AU). Consequently, spectra of giant planets at sep- 
arations analogous to those in the solar system (a > 5 
AU) can be observed only with high-contrast imaging. 
Since currently achievable sensitivities can resolve only 
systems that are hot, hence likely young, the planets are 
also likely formed in situ at the observed separations, and 
provide valuable probes for models of planet formation. 

Recent photometric observations of directly imaged 
planets are increasingly highlighting the need for a new 
class of atmospheric models to explain the data. Inter- 
pretations of current observations are based on tradi- 
tional models that were used to interpret atmospheric 
spectra of brown dwarfs (Burrows et al. 2001, 2006; 
Baraffe et al. 2003; Saumon & Marley, 2008). While 
giant planets at large orbital distances, and, hence, with 
low stellar irradiation, could have been expected to have 
had similar spectral characteristics to low temperature 
brown dwarfs, recent observations are proving otherwise 
(Marois et al. 2008 & 2010, Bowler et al. 2010, Currie 
et al. 2011). 

The HR 8799 system (Marois et al. 2008,2010) with 
four planetary companions at wide orbital separations is 
a case in point, epitomizing the difficulty in modeling di- 
rectly imaged planets. In the discovery paper, Marois et 
al. (2008) reported six-channel near-infrared photome- 
try of HR 8799 b, c, and d, and derived effective tem- 
peratures (T G ff), masses (M), and radii (R p ) from the 
observed luminosities and an assumed common age of 
the objects, using already published theoretical cooling 
tracks (Baraffe et al. 2003). On the other hand, deriving 
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the atmospheric parameters from spectral fits to their 
photometry, instead of cooling tracks, yielded unrcalisti- 
cally small radii (see Supplementary Online Material of 
Marois et al. 2008). This discrepancy was attributed to 
missing atmospheric physics, including inadequate cloud 
modeling, suggesting the need for new spectral models. 

Several follow-up observations of the HR 8799 system 
have also highlighted the need for new models. Bowler 
et al. (2010) reported spectroscopy of HR 8799b in the 
2.12 - 2.23 /mi range. While existing models were un- 
able to explain all their observations, the better fits were 
obtained with models with thick clouds. Furthermore, 
the T e ff and R p they derived from the best fitting spec- 
tral model were inconsistent with those derived based on 
evolutionary models. Hinz et al. (2010) reported pho- 
tometry at 3.3, 3.8, and 4.8 /xm, and suggested that mod- 
els with non-equilibrium chemistry provided better fits 
to their data, while noting the differences in fits based 
on different model assumptions. Based on VLT spec- 
troscopy in the 3.88-4.1 /jm range (Janson et al. 2010), 
similar conclusions were drawn for HR 8799c. 

Currie et al. (2011) further explored the HR 8799 plan- 
ets' atmospheres with data from Marois et al.(2008), new 
Subaru and VLT data, and re-reduced MMT data first 
presented in Hinz et al. (2010). Using the model formal- 
ism of Burrows et al. (2006), they found that standard 
L/T dwarf atmosphere models over-predicted the flux 
between 1 and 1.25 /xm. They found that "thick cloud" 
models, where the cloud density tracks that of the gas 
high in the atmosphere, reproduced the suppressed 1— 
1.25 micron flux and, thus, provided far improved fits 
to the data. However, they suggested that these mod- 
els may "overcorrect" for the weak cloud coverage in 
the standard models, as their preliminary 'partly cloudy' 
models fit the data better. They also acknowledged that 
their grid sampling was too coarse to precisely deter- 
mine best-fit surface gravities and effective temperatures 
from the data. Thus, further refinements of the "thick 
cloud" formalism with a new, expansive, set of atmo- 
sphere models are necessary to more accurately model 
data for known planet-mass objects and predict proper- 
ties of soon-to-be imaged planets. Our present work is 
motivated by this need. 

In this work, we embark upon an extensive modeling 
effort in order to present a new grid of atmospheric mod- 
els to aid in the characterization of directly imaged plan- 
ets like those of the HR 8799 system. We report models 
with different cloud types, both in physical extent and 
cloud composition, computed over a large grid in effective 
temperature, gravity, and metallicity. The models span 
a wide range of possible cloud cover, from clear atmo- 
spheres to those with clouds extending all the way to the 
top of the atmosphere, and encompass models that have 
been known to be consistent with field brown dwarfs. We 
follow the cloud parametrization developed in Burrows et 
al. (2006). Other studies in the past have also consid- 
ered alternate cloud models incorporating a wide range 
of geometrical and optical thicknesses (e.g. Marley et al. 
2002; Tsuji 2002, 2005; Helling et al. 2011). 

To demonstrate the accuracy of our models, we use the 
HR 8799 planetary system as a test case. We explore the 
model parameter space in search of the best fitting mod- 
els to the photometric observations of HR 8799b, c, and 
d, while also discussing in detail the degeneracies be- 



tween the various model parameters. As anticipated in 
recent literature (Marois et al. 2008, Bowler et al. 2010, 
Currie et al. 2011), the observations of the HR 8799 plan- 
ets indicate a potentially new class of objects with cloud 
properties unlike those of any field brown dwarfs known. 
We do not consider the recently discovered fourth com- 
panion HR 8799c (Marois et al. 2010), for which data are 
presently inadequate for a detailed atmospheric analysis. 

We report constraints on the atmospheric properties, 
masses, and ages of HR 8799b, c, and d. We use existing 
observations of near-infrared photometry for these ob- 
jects compiled in Currie et al. (2011). Our best-fit mod- 
els allow us to constrain T e s, log 10 (g), mass (M), and age 
(t) for each object. The range of estimated masses for 
HR 8799b, HR 8799c, and HR 8799d conservatively span 
2-12 Mj, 6-13 Mj, and 3-11 Mj, respectively. The best 
fits for HR 8799b, c, and d yield T cS of -850 K, ~1000 
K, and ~900 K, respectively, with log 10 (g) of 4.3, 4.2, 
and 3.8, respectively. The corresponding estimates of 
the masses are ~12, ~11, and ~6 Mj, respectively, and 
of the ages are ^150, ^65, and ~20 Myr, indicating that 
they are indeed roughly coeval, if the adopted age of the 
star is ~30-160 Myr. 

The observations of all the objects can be explained 
with models involving forstcrite clouds, with physically 
thick cloud cover, and 60-/im modal grain size. Models 
without clouds cannot explain the observations for any 
of the objects. The observations are not explained by 
canonical model spectra of L and T brown dwarfs (Bur- 
rows et al. 2006; Saumon & Marley 2008) at the same 
T G ff j indicating that the present objects might form a 
new class with geometrically thicker clouds than those 
inferred for L and T dwarfs, as also suggested in Cur- 
rie et al. (2011). The recent inference of clouds in the 
late T Dwarf Ross 458C (Burgasser et al. 2010) suggests 
that this object may be related. While the requirement 
of clouds at altitude is stringent, the cloud composition 
and particle sizes are less constrained by the data, as is 
the gas-phase metallicity. We find that models with Fc 
clouds, with l-/im size particles and 1% supcrsaturation, 
can also fit the infrared observations. Models including 
non-equilibrium chemistry do not obviate the require- 
ment for clouds. 

After establishing the accuracy of our models in con- 
straining the HR 8799 planet properties, we present a 
new suite of thick cloudy atmosphere models to facilitate 
characterization of future detections of similar objects. 
We report the physical properties and absolute photo- 
metric magnitudes, in seven near/mid infrared band- 
passes. Our model grid spans a wide range in planet 
properties, with T c g = 600 to 1700 K, and log 10 g = 
3.5 to 5.0, corresponding to masses between 2 and 40 
Mj, and ages between 4 Myr and 5 Gyr, and has cloud 
composition and structure similar to those of our best fit 
models for the HR 8799 planets. We explore the range in 
color-magnitude space occupied by our models and com- 
pare them with the locations of known field L/T brown 
dwarfs and other substellar objects. Our models predict 
that objects with thick clouds, like the HR 8799 planets, 
will occupy a distinguishable locus in color-magnitude 
space, diverging from the L/T dwarf sequence at the L/T 
transition, and extending to systematically fainter and 
redder regions in color-magnitude diagrams. 

In what follows, we first describe the atmospheric mod- 
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els in section [21 followed by a general exploration of the 
model parameter space in section [3j in order to elucidate 
the dependence of model spectra on the various param- 
eters. In section 21 we describe the data used in this 
work and the fitting procedure, and present model fits to 
observations of HR 8799b, c, and d. Motivated by the 
model solutions for the HR 8799 planets, we provide a 
new collection of models in section [5] to facilitate charac- 
terization of future detections of similar objects. Finally, 
in section [6l we summarize our conclusions. 

2. METHODS 

To generate atmospheres for this paper we employ the 
code COOLTLUSTY (Hubcny 1988; Hubeny & Lanz 
1995; Burrows et al. 2006,2008). We assume hydrostatic, 
radiative, and chemical equilibrium (generally) and ig- 
nore any irradiation by the central star. At the orbital 
distances of these objects this is an extremely good ap- 
proximation and breaks down only in the blue and ul- 
traviolet portions of the spectrum. We use the equa- 
tion of state of Saumon, Chabrier, & Van Horn (1995), 
the chemical code (from which we obtain the molecular 
abundances as a function of temperature and pressure) 
of Burrows & Sharp (1999) and Sharp & Burrows (2007), 
and the opacities compiled in Sharp & Burrows (2007) 
(see also Burrows et al. 2001). Our default spectral 
model employs 5000 frequency points from ^0.35 to 300 
microns. The major gas-phase species found in these 
theoretical atmospheres are H2, He, H2O, CO, CH4, N2, 
NH 3 , FeH, CrH, Na I, and K I. 

Refractory silicates condense out at temperatures be- 
low ^2300 K and above ^1700 K. The condensation of 
these refractories is automatically handled by the chem- 
ical code, which "rains out" the associated species when 
the condensation curve indicates, leaving a gas phase de- 
pleted of the constituent elements in stochiomctric ra- 
tios. The resulting elemental abundances of the depicted 
atmospheres are then used to determine the chemical- 
equilibrium molecular abundances. This procedure ac- 
counts for the settling of the refractories in a gravi- 
tational field and the depletion at altitude. However, 
clouds of these refractories form and assume shapes and 
thicknesses that reflect the meteorology of the atmo- 
sphere, in particular influenced by convection. The re- 
sulting cloud structures and grain properties (e.g., parti- 
cle sizes) are unknown and are the most difficult feature 
of brown dwarf and hot giant exoplanet atmosphere mod- 
eling. The thickness of the clouds likely scales in some 
way with the pressure scale height, but how is almost un- 
constrained by theory (Ackerman & Marley 2001; Helling 
et al. 2001,2004). 

Therefore, to model clouds we parametrize the spatial 
extent and particle size distributions. The approach we 
use for each is taken from Burrows et al. (2006) and 
Sudarsky et al. (2000), respectively. Because numerous 
refractory species condense out over a range of tempera- 
tures and the base of a cloud is likely to intersect the asso- 
ciated condensation curve ("Clausius-Clapcyron" equa- 
tion), we put the base of the clouds we study here at 
2300 K (approximately where the most refractory species 
condense). At the cloud base, it is assumed that the ele- 
mental abundances of silicon and magnesium (dominant 
elements in silicate refractories) or iron are solar, or when 
non-solar atmospheres arc considered that they scale ac- 



cordingly. Similarly, when we focus on iron clouds, the 
elemental abundance at the base of an iron cloud is as- 
sumed to be continuous across the cloud/gas interface. 
The upward extent in pressure from this base is then 
parametrized in terms of some multiple or fraction of 
the local pressure scale height. The baseline distribution 
of the cloud particle density is assumed to follow the 
gas-phase pressure profile and we multiply this baseline 
reference distribution by a shape function, /(P), to de- 
termine the actual model cloud distribution in pressure 
space (hence, real space). /(P) is never greater than one 
and can be made to cut off sharply at the top and bottom 
of the cloud. Specifically, we define /(P) as: 

r (p/p u ) Su , p < p u 

f{P) = < /cloud, Pu<P<Pd (1) 

{(P/P d )- Sd , p>p d , 

where P d is the pressure at the cloud base and P u < P d . 
In particular, P u is the pressure at the intersection of 
the atmospheric T/P profile and the condensation curve 
of the least refractory condensates, assumed here to be 
either forsteritc or iron. The indices s u and s d define the 
rapidity with which the clouds are cut off on their upper 
and lower boundaries. The larger they are, the sharper 
the cutoff. When / c i ou d = 1 the cloud has a flat portion 
in its middle between P u and Pj. The parameters s u , s d , 
and /cloud define the cloud spatial structure. 
For this study, we consider four model parameter sets: 

Model E: s u = 6, s d = 10, / cloud = 1 
Model A: s u = 0, s d = 10, / cloud = 1 
Model AE: s u = 1, s d = 10, / cloud = 1 
Model AEE: s u = 2, s d = 10, / c i oud = 1. 

Since s d = 10, all these models have a sharp cutoff at 
high pressures (at their base). Also, all these cloud mod- 
els have a "flat" profile in the middle that follows the 
pressure profile itself. The latter is roughly exponential 
with the pressure scale height. It is in the extent of the 
cloud at altitude that these models differ. Model E cuts 
off rapidly at altitude. Model A extends all the way to 
the top of the atmosphere, with a particle number density 
distribution that follows the gas. Model AE is interme- 
diate between Models E and A, decaying at the cloud 
tops with a scale height that is equal to twice that of the 
gas. Nevertheless, it still extends beyond Model E and 
puts particles higher up at lower pressures. Model AEE 
is similar to Model AE, but extends to altitude with a 
scale height of a third of a pressure scale height. Model 
E was the model settled upon in Burrows et al. (2006) 
as the best fitting model for L and T dwarfs, though it 
still did not capture the rapidity of the L-to-T transition. 
To model the HR8799 planets, we are forced to extend 
the cloud tops to lower pressures, thickening the clouds 
to a specified degree. We have thus introduced the 'AE' 
cloud models in this work. 

We refer to the AE-type clouds as 'thick clouds' in this 
work, implying that they are thicker in physical extent 
than the E-type clouds that are known to fit observa- 
tions of L dwarfs. For the same grain size distribution, 
the physical thickness naturally also enhances the optical 
thickness of the cloud. Our goal in the present work is to 
constrain the combination of physical extent and cloud 
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grain properties of clouds that best explain the observa- 
tions of the HR 8799 planets. Recent studies have advo- 
cated the need for such thick cloudy models to explain 
the HR 8799 observations (Bowler ct al. 2010, Currie et 
al. 2011). Several theoretical studies in the past have 
also considered cloud properties comparable to our thick 
cloud models (e.g. Chabrier et al. 2000; Allard et al. 
2001; Tsuji 2002 & 2005; Burrows et al. 2006). 

For the size distribution of cloud particles, we use the 
Deirmendjian distribution of cumulus water clouds on 
Earth (Sudarsky et al. 2000): 



dN 
da 



cxp 



-6 



(3) 



The modal particle size is ao and is a free parameter in 
our model. In this paper, we assume spherical particles 
of either forsterite (Mg2Si04) or iron, and use Mie scat- 
tering theory with a table of complex indices of refrac- 
tion as a function of wavelength for each (Sudarsky et al. 
2000). In Burrows et al. (2006), the fiducial best-fitting 
model had not only a Model E structure, but forsterite 
particles with a modal size of 100 microns. Since iron un- 
dergoes homogeneous condensation (i.e. vapor and liq- 
uid/solid phases composed of the same chemical species), 
we assume a supersaturation factor of 1% to represent 
the fraction of iron in the condensed phase (Fegley & 
Lodders, 1996; Ackerman & Marley 2001; Marley et al. 
2002; Burrows et al. 2006). As mentioned above, the 
iron elemental abundance at the base of an iron cloud is 
assumed to be continuous across the cloud/gas interface. 
Any elemental discontinuities would be smoothed out by 
upward gaseous diffusion and overshoot from the lower, 
hotter, regions of the atmosphere. 

The opacity spectrum of forsterite clouds is simi- 
lar to that of enstatite clouds in mid-IR wavelengths 
(A > 5/im). In the near-IR, the forsterite clouds have 
higher opacity in the 1—3 [im range compared to en- 
statite clouds. On the other hand, the opacity spectrum 
of Fe clouds is significantly different from both enstatite 
and forsterite clouds. The Fe clouds offer substantially 
higher grain opacities in the near-infrared, with a steep 
decline into the mid-infrared. We do not consider in de- 
tail the dependence of cloud formation processes on the 
gas density and gravity, although cloud formation is fa- 
vored at higher densities, and the microscopic nature of 
clouds in the different gravity regimes could be different. 

3. DEPENDENCE OF PREDICTED SPECTRA ON MODEL 
PARAMETERS 

Before describing detailed model fits to observations of 
the HR 8799 planets, we present a general exploration of 
the parameter space to understand some generic features 
of model spectra and their dependence on model param- 
eters. The parameters in our model are: T e gf, log 10 (g), 
metallicity, cloud composition, modal particle size, and 
the physical extent of cloud cover. We also explore mod- 
els of non-equilibrium chemistry, in which case the eddy 
diffusion coefficient, K zz , (Saumon et al. 2003; Hubeny 
& Burrows, 2007) is an additional parameter. 

3.1. Cloud-free models 

We first consider model atmospheres without clouds. 
Figure [TJ shows cloud-free models as a function of the 



various model parameters, along with observations of 
HR 8799b, shown for reference. Here, we focus on a 
nominal spectrum with T e g = 1000 K and log 10 g = 
4.5, the red model in Fig. [1] The model nomencla- 
ture specifies the atmospheric parameters. For example, 
'T1000.1s.g4.5' refers to a cloud-free (T) model with 
T e ff = 1000 K, log 10 (g) = 4.5, and lxsolar metallicity. 
The major features in the spectrum can be understood in 
terms of the dominant opacity sources. The high fluxes 
in the J, H, and K bands (within ~1 — 2.4 /j,m) are due 
to windows in water vapor opacity, and arc a defining 
characteristic of cloud-free models; the gaps between the 
high fluxes correspond to absorption bands of water va- 
por. The deep absorption trough between 2.2 and 3.5 
/zm is due to water and methane absorption. Water ab- 
sorption also causes deep features between 5 and 7 /xm, 
while the peak between 4 and 5 /xm is due to the lack 
of significant absorption features due to any of the dom- 
inant species. In principle, CO has strong absorption 
features between 4 and 5 /xm. However, for T e g = 1000 
K, and assuming equilibrium chemistry, methane is the 
dominant carbon carrier, as opposed to CO. 

The effects of T ff , g, metallicity, and non-equilibrium 
chemistry on model spectra of cloudless atmospheres are 
shown in Fig. [TJ An increase in T e ff causes an increase 
in emergent flux across all wavelengths. However, the 
relative increase in flux varies with wavelength, depend- 
ing on the molecular absorption. In regions of low water 
opacity, such as in J, H, and K, and between 4 and 5 /xm, 
the emergent flux is highly sensitive to temperature. On 
the other hand, the emergent flux is less sensitive to tem- 
perature at wavelengths where the spectrum is saturated 
with water absorption, as shown in the top-left panel of 
Fig. [TJ It is also apparent that none of the cloudless mod- 
els can explain the observations of HR 8799b, shown for 
reference, at wavelengths below 2.2 /im. 

The dependence of spectra on gravity is shown in the 
top-right panel of Fig. [TJ A decrease in gravity causes 
a general increase in flux across all wavelengths, except 
in regions saturated by strong water vapor absorption. 
The dependence of the spectra on gravity arises from 
two factors. For a given T e g, a lower gravity causes an 
increase in radius (see Eq. 5 of Burrows et al. 2001), 
leading to a corresponding increase in the emergent flux. 
In addition, a higher gravity causes a higher mean density 
in the atmosphere, leading to higher opacities, especially 
in the H2-H2 collision-induced opacity, which increases 
with density. The higher opacities then cause enhanced 
flux suppression in all but the 3.3-/xm channel which is 
saturated with methane and water opacity. 

At the level of present observations, and in the context 
of our models, the spectra are not particularly sensitive 
to changes in metallicity or to non-equilibrium chem- 
istry due to eddy mixing. As shown in the bottom-left 
panel of Fig. [TJ variation in metallicity between solar and 
3 x solar values, does not significantly alter the spectra in 
the observed bands. Finally, as shown in the bottom- 
right panel of Fig. [TJ non-equilibrium chemistry induces 
nominal changes in the model spectra, albeit too weak 
to be determined with current observations. For mod- 
els with non-equilibrium chemistry, the eddy diffusion 
coefficient (K zz ) causes transport of CO from the hot- 
ter lower atmosphere to the upper, observable regions of 
the atmosphere (Saumon et al. 2003; Hubeny & Bur- 
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rows 2007). Consequently, higher K zz and faster diffu- 
sion leads to more CO absorption, which is apparent in 
the pronounced absorption features in the 4-5 /im spec- 
tral region. 

3.2. Cloudy models 

Model spectra of atmospheres with clouds are shown 
in Fig. [2] Four types of cloudy models (E, A, AEE, and 
AE) are shown, along with observations of HR 8799b, 
shown for reference. As described in section [21 the A 
models have cloud cover extending all the way to the top 
of the atmosphere, whereas the E models cut off sharply 
at altitude. The AE and AEE models are intermediate 
between A and E models, with the AEE models lean- 
ing more towards the E models for the same T e s and 
g. All the models have equilibrium chemistry with so- 
lar abundances. The model nomenclature encapsulates 
the cloud type and atmospheric parameters. For exam- 
ple, 'E60.1000.1s.g4.0' implies 'E' type cloud cover with 
60-/xm modal particle size, T e g = 1000 K, 1 x solar metal- 
licity, and log 10 g = 4.0, where g is in cm s~ 2 . 

Model spectra of atmospheres with 'E' type clouds are 
shown in the top-left panel of Fig. [2] E-type models have 
been shown to provide good hts to observations of L and 
T dwarfs (Burrows et al. 2006). Because E-type models 
have minimal cloud cover among the models considered 
here, their spectral features bear some resemblance to 
those of cloud-free models shown in Fig. [1] The large 
peaks in the z, J, H, and K bands show regions of low 
water opacity, as also seen in the cloud-free models. But, 
the peaks can be substantially lower, and troughs higher, 
compared with those of cloud-free models, at high T e g 
and g (Burrows et al. 2006). In the present context, 
however, the model fluxes in the z, J, and H bands are 
still too high by a factor of 2 - 5 than the observations of 
the HR 8799 objects in these bandpasses, as suggested 
by data for HR 8799b in Fig. [3 

Clouds of thicker physical extent than the E models 
cause greater suppression of the emergent flux in the z, 
J, and H bands. This situation is seen in models with A- 
type clouds, as shown in the top-right panel of Fig. [21 
While A models can, in principle, yield low observed 
fluxes in the z, J, and H bands, they predict much higher 
fluxes than E models at longer wavelengths, in the L 
and M bands, for the same T e g. This is caused by the 
spectral redistribution of energy from the shorter wave- 
lengths, where cloud opacity suppresses the emergent 
flux, to longer wavelengths in the spectrum. However, 
as shown in the top-right panel of Fig. [21 A models can- 
not fit all the observations of HR 8799b simultaneously. 

The above discussion suggests that cloudy models with 
physical extent intermediate between A and E models 
might be able to explain the data better. The bottom- 
left and bottom-right panels of Fig. [2] show such mod- 
els, i.e., the AEE and AE models described in section [21 
While AEE models, those that are closer to the E mod- 
els, are still inadequate to explain the observations, the 
AE models provide the best match to the data. The 
thicker physical extent of the clouds in the A and AE 
models over those of the E model are apparent in the 
total optical depths contributed by each to the observ- 
able atmosphere. For example, considering the models 
in Fig. [2] with T cff = 1000K and log 10 g = 4.0, the total 
Rosseland optical depths at pressures below 1 bar for the 



A, AE and E models are 12.6, 5.0, and 1.2, respectively. 
Considering pressures below 0.3 bar, the total Rosseland 
optical depths for the same A, AE, and E models are 4.4, 
1.1, and 0.6, respectively. In what follows, we present 
the best fits to the observations of the three objects 
HR 8799b, c, and d, based on AE models. Figure [3] shows 
temperature-pressure (T-P) profiles for four representa- 
tive cloudy models. A-type forsterite cloud models which 
have the thickest (physically) cloud cover require hotter 
temperature profiles, for the same T c ff = 1000 K and 
l°gio5 = 4-0, compared to AE-type and E-type models. 
On the other hand, AE-type forsterite clouds and A-type 
iron clouds, with 1% supersaturation, yield similar T-P 
profiles. 

4. APPLICATION TO THE HR 8799 PLANETS 

In this section, we present model spectra that fit the 
observations of HR 8799b, c, and d, along with estimates 
of masses and ages. For each object, we have extensively 
explored the parameter space of models with and without 
clouds in search of the best-fitting models. As described 
in section [3] for the case of HR 8799b, we find that cloud- 
free models cannot explain the data for any of the three 
objects, and that cloudy models with AE-type clouds 
provide the best match to the data. Below, we present 
the constraints on the T e g, log 10 <?, mass (M), and age 
(t) for each object, along with the best fitting models. 

4.1. Observations and Fitting procedure 

Several recent studies of the HR 8799 system have re- 
ported broadband photometry and colors of the compan- 
ions. For the spectral fits in the present study, we use 
broadband photometry at nine wavelengths reported in 
the literature for HR 8799 b, c, and d. The observa- 
tions we use are described and summarized in Currie et 
al. (2011). The nine filters and their central wavelengths 
are: z (1.03 /imfl J (1.248 /urn), H (1.633 /im), CH4 S 
(1.592 /im), CH4 L (1.681 /im), K s (2.146 /im), [3.3] (3.3 
/im), L' (3.776 /im), and M (4.75 /im). The photometry 
in the J, H, CH4s, CH4l, and K s bands for all the three 
objects is taken from Marois et al. (2008). The obser- 
vations in the [3.3], L', and M bands are adopted from 
Currie et al. (2011). The observations in the M-band 
for all the three objects are 3-er upper limits, as are the 
z-band observations of HR 8799c and d, and the [3.3]- 
band observation of HR 8799d. We do not consider the 
recently discovered fourth companion HR 8799e (Marois 
et al. 2010; Currie et al. 2011) as data for this planet 
are currently restricted to Ks and L' band and thus are 
inadequate to place meaningful constraints on its atmo- 
spheric parameters. 

We determine the best-fit models to observations of 
each object using a grid-based approach. We computed 
models over a large grid in T c ff and log 10 g for differ- 
ent cloud types, particle sizes, and metallicities. After 
a global exploration of the grid, we found that AE type 
cloudy models provided the best match to the data. Con- 
sequently, we report error contours in a two-dimensional 
(T e ff - log 10 g) grid of AE-type forsterite cloud models. 

3 The Subaru/IRCS z band filter used by Currie et al. (2011) is 
very similar to the better known MKO "Y" band filter. The IRCS 
z filter and Y filter should not be confused with the Sloan z' filter, 
which covers shorter wavelengths. 
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We nominally fix the modal particle size at 60 /im and 
the metallicity at solar, but we found that forstcrite par- 
ticle sizes of 30 - 100 /im provided almost equally good 
fits. We evaluate the model fits using the x 2 statistic, 
given by: 

N obs / , f \ 2 

2 X / Ji,model Ji,obs \ / ,\ 



where, fi.modei is the channel-integrated emergent flux of 
the model in each channel, and fi i0 b s and <7i. 6s are the 
observed flux and the l-a uncertainty in that channel, 
respectively. N t, s is the number of observations. Where 
only 3<7 upper limits are available, we compute their con- 
tribution to the x 2 using the following approach. When 
the model flux is less than the upper limit, we do not 
penalize the \ 2 - When the model flux exceeds the upper 
limit, the contribution to the \ 2 is {^[f model - f+}/f+) 2 , 
where /+ is the observed 3cr upper limit. 

We derive the confidence limits on the model fits based 
on the x 2 distribution and the number of degrees of free- 
dom (Ndof)- When N b s = 9 observations (Ndof = 8) 
are used in computing the x 2 ■, the 3er, 4<7, 5<7, 6er, and 
7a confidence levels correspond to x 2 of 23.6, 33.0, 44.0, 
56.8, and 71.1, respectively. For N obs = 8, i.e. N do f = 7, 
the corresponding values of x 2 are 21.8, 31.0, 41.8, 54.4, 
and 68.5, respectively. 

4.2. HR 8799b 

Based on a detailed exploration of cloudy models, we 
find that the observations of HR 8799b are best explained 
by models with forsterite clouds with physically thick 
AE-type cloud cover, and a nominal modal particle size 
of 60 /im. As described in section[2j AE-type cloud mod- 
els have cloud cover that is intermediate between the A 
and E type cloud models. As shown in the top-left panel 
of Fig. SI the observations can be explained with an AE- 
type forsterite cloud model with T e g = 850 K and log 10 g 
= 4.3. The observations can also be explained by models 
with thick A- type Fe clouds with 1-fim modal particle 
size and 1% supersaturation, implying multiple possibil- 
ities for the cloud composition and particle size. In the 
remaining three panels of Fig. 21 we show the depen- 
dence of our best-fit model on T g, log 10 g, and cloud 
properties. An increase in T e g or a decrease in g over the 
best-fit model leads to systematically higher fluxes across 
the spectrum, as described in section [3] The bottom- 
right panel shows the dependence on cloud-cover type. 
A- type clouds, with complete cloud cover up to the top 
of the atmosphere, cause excess suppression of flux in 
the 1-2 /im range. On the other hand, E-typc clouds 
that have cloud cover truncated at height are insufficient 
to suppress the 1-2 /im flux to the observed levels. The 
degeneracy in cloud composition, between forsterite and 
Fe clouds, might be resolved by future observations in 
the visible and mid-infrared. Iron clouds have strong 
absorption in the visible, whereas silicates have strong 
spectral features at ^10 /im and ^20 /im which might 
be observable in high resolution spectra (Cushing et al. 
2006; Helling et al. 2006; Henning 2010). 

The space of models explored in order to fit the obser- 
vations of HR 8799b, along with distribution of the x 2 in 
this space, are shown in Figs. [5] and [6] Each model in the 



large grid explored is shown color-coded by the degree of 
fit to the data, denoted by the confidence level corre- 
sponding to the x 2 ■ The purple open-diamonds show 
fits to within the 3<7 confidence level, and the red, blue, 
green, black, and grey, filled-circles indicate fits within 
4cr, 5cr, 6cr, 7cr, and > 7<7, confidence levels, respectively. 
In Fig. [51 the x 2 is computed by taking into account 
all the nine photometric observations described in sec- 
tion 14.11 In this case, we obtain a high x 2 f° r all the 
models, with only a few models lying within the 3<7 level. 
We find that a large part of the x 2 1S contributed solely 
by the 3.3-/zm data point, as is also evident from Fig. |4l 
This discrepancy can in principle be explained by a lower 
methane abundance than the equilibrium abundances (at 
solar metallicity) considered in our cloudy models. With- 
out the contribution of the 3.3-/xm point to the x 2 we 
obtain much better fits to all the remaining data, while 
still preserving the best-fit model parameters, as shown 
in Fig. [6] Instead of generating an entirely new grid of 
models with less methane in order to fit the single 3.3-/zm 
point, we assess the model fits based on the remaining 
eight observations. We consider the latter procedure as 
a more accurate means of assessing the model fits to the 
data and quote the ranges of fit parameters accordingly. 
Consequently, we have N b s = 8 observations, yielding 
Ndof = 7 degrees of freedom. 

The ranges of the physical parameters of HR 8799b at 
different levels of confidence are shown in Fig. [5] The 
constraints at 3er confidence and the parameters for the 
best-fit model are shown in Table Q] and Table [H The 
solutions at 3a confidence yield T c s = 750 - 850 K, log 10 g 
= 3.5 - 4.3, M = 2 - 12 Mj, and t = 10 - 150 Myr. 
The best-fitting model, reported in Table [lj has T e g = 
850 K and log 10 g = 4.3. From the T e g and log 10 .g, we 
derive a mass of ~12 Mj and an age of ^150 Myr, using 
analytic fits to evolutionary models in the present grid 
(e.g. Burrows et al. 2001). 

4.3. HR 8799c 

Our model fits to observations of HR 8799c require 
cloud properties that are similar to those of HR 8799b, 
but with different T c ff and log 10 g. Figure [7] shows model 
spectra of HR 8799c and their dependence on the model 
parameters. An AE-type forsterite cloud model with T ff 
= 1000 K and log 10 g = 4.2 provides a good match to the 
observations. The observations can also be explained 
with Fe clouds with l-/tm modal particle size and 1% 
supersaturation, but with similar T c ff and log 10 g of 1000 
K and 4.0, respectively. The dependence of the best-fit 
model spectrum on T e g, log 10 <7, and cloud cover follow 
the general trends described for HR 8799b above. 

Figures [8] and [9] show the grid of models explored for 
HR 8799c color-coded with the level of fit to the data. 
We find that considering all the nine observations yields 
large values of x 2 , with no model in the grid fitting to 
within 3cr confidence. As with HR 8799b, we find that 
the observed 3.3-/im data point contributes the most to 
the x 2 ■ The flux observed in the 3.3-/im band is higher 
than the model predictions in this bandpass by over 2- 
(7, and suggests methane depletion in HR 8799c, as also 
suggested for HR 8799b. Consequently, following the ap- 
proach discussed in section 14.21 we assess the model fits 
without the 3.3-/im point, as shown in Fig. [9l and quote 
the fit parameters accordingly. 
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The constraints on the model parameters at different 
levels of confidence are shown in Fig. [9] and Table [2] At 
3cr confidence, the solutions for HR 8799c yield T e g = 
925 - 1025 K, log lo5 = 3.8 - 4.3, M = 6 - 13 Mj, and t 
= 20 - 100 Myr. A representative best-fit model has T e g 
= 1000 K, log 10 g = 4.2, M = 11 Mj, and t = 65 Myr, as 
shown in Table[T] In principle, an equally good fit is also 
obtained for T c g = 1050 K, and log 10 g = 4.5, as shown in 
Fig. [6] in the purple diamond bounded by a black square, 
yielding a mass of 17 Mj. But we rule out the latter 
solution since the corresponding age of the object would 
be 180 Myr, which is larger than the maximum age of 
160 Myr estimated for the host star. 

4.4. HR 81996, 

The observations of HR 8799d can be explained by a 
broad set of model parameters. Cloudless models are 
still inadequate to explain the data, as with HR 8799b 
and HR 8799c. However, a wider range of cloudy mod- 
els fit the observations. As shown in the top-left panel 
of Fig. [T0l forsterite clouds of AE type and A type can 
explain the observations, with different T g and log 10 g. 
The best-fit model with AE-type cloud cover has a T e g 
= 900 K and log 10 g = 3.8 as shown in Table [T] On the 
other hand, a model with A-type cloud cover with T e g 
and log 10 g of 1000 K and 4.0, respectively, also provide 
a reasonable match to the data which comprises of three 
upper-limits and are, hence, less constraining. The de- 
generacy in cloud type leads to multiple estimates of the 
best-fit mass and age of HR 8799d. The best-fit model 
with AE-type clouds yields M ~ 6 Mj and t ~ 20 Myr, 
whereas, that with A-type clouds yields M ~ 8 Mj and 
t ~ 30 Myr. With this level of degeneracy, the obser- 
vations are insensitive to slight changes in metallicity or 
to the presence of non-equilibrium chemistry. The obser- 
vations can also be explained by models with Fe clouds, 
1-^m modal particle size, 1% supersaturation, and T e g 
and log 10 .g of 1000 K and 4.0, respectively. The depen- 
dence of the best-fit models on the model parameters is 
shown in Fig. [TO] 

The range of fitting AE models for HR 8799d is shown 
in Fig. [TT]at different levels of fit, and the 3er constraints 
are reported in Table [5] At 3<r confidence, the ranges in 
the different physical parameters are: T g = 850 - 975 
K, log lo5 = 3.5 - 4.2, M = 3 - 11 Mj, and t = 10 - 70 
Myr. We note that, unlike HR 8799b and HR 8799c, the 
presence or absence of the 3.3-/xm point does not signif- 
icantly influence the \ 2 for model fits to the HR 8799d 
data. This is because for HR 8799d only an upper limit 
was observed at 3.3 fim, which is higher than the flux 
predicted in this channel by all the fitting models. 

5. NEW MODELS FOR THICK CLOUDY ATMOSPHERES: 
PREDICTIONS FOR 1-5 /iM SPECTRA 

Our model fits to observations of the HR 8799 plan- 
ets hint at a new regime in substellar atmospheres. The 
clouds required to explain the observations of HR 8799b, 
c, and d, have much thicker physical extent compared 
to those traditionally required to explain the observa- 
tions of L and T brown dwarfs at similar effective tem- 
peratures. The high-altitude clouds we infer in the 
HR 8799 planets might be a result of their lower grav- 
ities (log 10 g ~ 3.5 — 4), compared to most measured L 
and T dwarfs (Burrows ct al. 2006), which can lead to 



slower gravitational settling of condensate particles (e.g. 
Helling 2011). Additionally, the thick clouds might also 
be indicators of high metallicity and significant vertical 
mixing, which might be constrained by future data. 

The observed fluxes of HR 8799b, c, and d in the near- 
infrared wavelengths (in z, J, H, and Ks bands) are two to 
five times lower compared with model spectra of T dwarfs 
at similar T e g. Major direct imaging programs currently 
underway - e.g. IDPS and Gemini/NICI (Marois et al. 
2011, in progress; Liu et al. 2010) - and upcoming pro- 
grams like GPI and SPHERE (Macintosh et al. 2008; 
Bcuzit et al. 2008) will likely image many new planets 
with a wide range of ages and masses, which translate 
into a wide range in T e g and g (e.g. Burrows et al. 1997, 
2001). Because the mid-IR portion of a planet's SED 
is strongly affected by T c g and g, and the presence of 
clouds, a wide spread in photometric colors and magni- 
tudes can be expected from future detections. 

We are, therefore, motivated to provide a suite of 
model spectra in order to aid in the characterization of 
objects like the HR8799 planets likely to be detected in 
the near future. We tentatively refer to such objects as 
those bearing 'thick-cloud' atmospheres, and give them 
the AE-type cloud properties required to explain the HR 
8799 objects. Such models might also help in the charac- 
terization of objects like the cloudy T Dwarf Ross 458C 
(Burgasscr ct al. 2010). We present in Table [3] mod- 
els of thick-cloud type atmospheres over a representative 
grid in T e g and log 10 g. For each model in the grid, we 
provide mass (M), age (t), and absolute magnitudes in 
seven photometric bandpasses in the infrared (z', Y, J, 
H, Ks, L', and M). The model grid spans T eB = 600 to 
1700 K, and log 10 <? = 3.5 to 5.0, where g is in cm s~ 2 . 
The corresponding masses range between 2 and 40 Mj, 
and ages range between 4 Myr and 5 Gyr. We assume 
the atmospheres to be in chemical equilibrium with solar 
abundances. 

5.1. Correlations between g/r c e and Mass/ Age 

Low-mass objects of the thick-cloud type occupy a dis- 
tinct region in atmospheric phase space. As shown in 
Table [3J objects with masses below 10 Mj are charac- 
terized by T off < 1100 K, log 10 5 < 4.5. On the other 
hand, higher mass objects (M > 15 Mj) are typically 
associated with higher specific gravities and/or higher 
T c g. This correlation can be understood by considering 
the dependence of mass on T e g and g (Burrows et al. 
2001). Mass has an approximate power- law dependence 
on both g and T e g, the dependence on g being steeper. 
For example, from Burrows et al. (2001), 

M ~ 35 Mj ( s 9 , \ ( eff T ) . (5) 
J V10 5 cm/s 2 / V1000 K/ V ' 

The age, on the other hand, has a power-law dependence 
on g, and a steeper inverse power-law dependence on T c ff . 

* i r ( 9 y-y ioooKy .s 

Consequently, young planetary objects are likely to have 
high T c ff and low surface gravities. These scaling rela- 
tions justify our choice of focusing on the subset of our 
model space with log(g) < 5. 
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The effect of T e ff and g on the emergent spectra are 
apparent from the infrared magnitudes. As shown in 
Tabic [3] and in Fig. |4l higher T e g and lower g result in 
systematically higher fluxes in all the bandpasses. How- 
ever, as shown in Fig. 21 the influence is highest in the 
Ks and L' bands, which are not saturated by water vapor 
absorption. 

5.2. Comparisons with Colors and Magnitudes of 
L/T-Type Substellar Objects 

To compare the predicted near-IR properties of objects 
with thick clouds with those with more standard cloud 
prescriptions, we generally follow the presentation struc- 
ture in Currie et al. (2011). We include color- magnitude 
loci for Case E standard cloud deck models from Burrows 
et al. (2006) for log 10 <? = 4 and 4.5 and a solar abun- 
dance of metals. To bracket the range in surface gravi- 
ties for planet-mass objects from our AE model grid, we 
overplot color-magnitude loci for the log 10 <? = 3.5 and 
4.5 cases. 

For empirical comparisons, we primarily use the L/T 
dwarf sample compiled by Leggett et al. (2010). To this 
sample, we add photometry for HR 8799 presented in 
Currie et al. (2011) and photometry for other planet- 
mass objects and very low-mass brown dwarfs with T e ff 
= 800-1800K: 2M 1207b (-5 Mj), 1RXJ1609. 1-210524 
(-9 M, 7 ), AB Pic (-13.5 Mj), and HD 203030b (-25 
Mj) (Chauvin et al. 2004, 2005; Mohanty et al. 2007; 
Lafreniere et al. 2008, 2010; Ireland et al. 2010; Metchev 
and Hillenbrand 2006). 

Figure [T^] presents our color-magnitude diagrams - 
H/H-K s , J/J-K s , J/J-H, and H/Y-H. In all four panels, 
the locus for the thick cloud models follows an extension 
of the L dwarf sequence to systematically fainter absolute 
magnitudes and redder colors than the standard cloud 
deck models. Both the cloud deck and thick cloud layer 
models loci pass through/near the color-magnitude space 
for L dwarfs and the hottest, youngest planet-mass ob- 
jects (1RXJ1609. 1-210524 and AB Pic), though the lat- 
ter's slightly redder colors are somewhat better matched 
by the thick cloud models. However, the standard mod- 
els provide far better matches to the colors and magni- 
tudes of T dwarfs. Conversely, the thick cloud layer mod- 
els reproduce the photometry of some objects near the 
L/T dwarf transition, but provide far better matches to 
the photometry for faint planet-mass objects/low-mass 
brown dwarfs, such as the HR 8799 planets, HD 203030b, 
and 2M 1207b. 

The differences between the models are most strik- 
ing at temperatures less than —1300 K, where standard 
cloud deck models substantially differ by predicting that 
progressively fainter/cooler substellar objects have bluer 
near-IR colors. For instance, for an object with log 10 g = 
4.5, T cff = 1200 K (w 17 Mj) the predicted position in 
J/J-K s is ps [14, 1] in the standard cloud deck case, but 
[14.3, 2.2] for the Model AE case. At lower temperatures 
(?cff ^ 1000 K) corresponding to objects with masses 
approaching or below the deuterium-burning limit, the 
differences are even more pronounced. 

5.3. Prediction of a Separate "Planet Locus" 
Identifiable from Upcoming Imaging Surveys 

Our models predict that objects with thick clouds and 
temperatures less than ps 1300 K will occupy far different 



regions of near-IR color-magnitude diagrams than those 
with clouds more characteristic of field brown dwarfs. 
As inferred from Table [3l objects with T c h < 1300 K and 
log 10 g = 3.5 and 4.5 will have masses less than —4 Mj 
and —18 Mj, respectively. The ages corresponding to 
objects with these ranges in surface gravity and temper- 
ature are >1 Myr and >117 Myr, respectively. 

Upcoming imaging surveys like GPI are sensitive to 
detecting 1-10 Mj around numerous stars with ages — 
10-100 Myr and, thus, will probe a range of temperatures 
and surface gravities where cloud deck and thick cloud 
models predict very different near-IR SEDs. Therefore, 
our analysis predicts that the near-IR SEDs of many 
soon-to-be imaged planets will not track the L/T dwarf 
sequence. Rather, like the HR 8799 planets and 2M 
1207b, they will track a distinguishable "planet locus" 
following an extension of the L dwarf sequence, diverg- 
ing from the L/T dwarf sequence for field brown dwarfs 
at the L/T transition, and extending to systematically 
fainter and redder color-magnitude positions. 

6. CONCLUSIONS 

We have generated new models to help characterize the 
atmospheres of the directly-imaged planets HR 8799b, c, 
and d, based on observations of broadband infrared pho- 
tometry. Recent observational efforts have led to good 
data on the HR 8799 system. However, explaining the 
observations with existing atmospheric models had thus 
far proven challenging. In the discovery paper, Marois 
et al. (2008) found that spectral fits to the photomet- 
ric observations yielded very high T e ff, or unrealistically 
small radii, for the objects, inconsistent with the values 
obtained from fits to cooling tracks. Subsequent follow- 
up observations (Bowler et al. 2010, Hinz et al. 2010, 
Janson et al. 2010, Currie et al. 2011) have all reported 
the inadequacy of existing models in simultaneously rec- 
onciling all the available observations. 

The new models of thick-cloud atmospheres we have 
generated in this work provide good fits to the observa- 
tions of HR 8799b, c, and d. Several recent studies on 
the atmospheres of these objects had all suggested the 
requirement of thick clouds to explain the observations 
(Marois et al. 2008, Bowler et al. 2010, Currie et al. 
2011). However, the lack of detailed models spanning 
different cloud parameters (physical thicknesses, cloud 
compositions, and grain size) had been an impediment to 
determining the nature of clouds required to fit the data. 
In a companion paper to this study, Currie et al. (2011) 
combined their new data with data from Marois et al. 
(2008) and a re-reduction of MMT data from Hinz et al. 
(2010) to identify models with intermediate cloud cover 
that might provide better fits to the observations. Their 
conclusions were based on our preliminary 'patchy' cloud 
models, derived using the modeling procedures found in 
Burrows et al. (2006). The present paper is meant to be 
a more detailed follow-on study to constrain the atmo- 
spheric parameters of the HR 8799 planets and to provide 
new models for future use. 

Our primary finding is that a thicker cloud cover than 
that required to explain observations of L and T brown 
dwarfs (e.g., the 'E' model of Burrows et al. 2006) is es- 
sential to explain the present observations of HR 8799b, 
c and d. Our best-fit solutions follow from an explo- 
ration of cloudy models with different physical extents, 
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cloud compositions, and modal particle sizes. In terms 
of the physical extent, we consider four cloud types: E, 
A, AE, and AEE. While E-type clouds truncate at alti- 
tude, A-type clouds extend all the way to the top of the 
atmosphere. The best-fit models for HR 8799b and c re- 
quire clouds intermediate between the A and E models, 
denoted as AE models in the present work; HR 8799d 
allows solutions with both AE and A models. Models 
with AE type forsterite clouds with a modal particle size 
of 60 /im were able to explain the observations of all the 
objects. However, A-type clouds made of pure Fe with 
l-/itm particle size and 1% super saturation could also ex- 
plain the observations reasonably well, implying multiple 
solutions for the cloud composition and particle sizes that 
can yield the required opacities. Despite the similar fits 
obtained with the different cloud compositions, we note 
that the opacity spectra of the forsterite and Fe clouds 
are significantly different. Future observations in the vis- 
ible and mid-IR might be able to distinguish between the 
two compositions. While forsterite clouds would be ex- 
pected to show silicate features near ~10 /im and ^20 
/im, Fe clouds have notable opacity in the visible. 

We derive constraints on T e g, g, mass (M), and age (t) 
for HR 8799b, c, and d, based on spectral fits to previ- 
ously published broadband photometric data in the near- 
infrared. At 3-cr confidence, the constraints for HR 8799b 
are T cff = 750 - 850 K, log 10 g = 3.5 - 4.3, M = 2- 12 
Mj, and t = 10 - 150 Myr. The parameters of the best-fit 
solution are T off = 850 K, log 10 g = 4.3, M ~ 12 Mj, and 
t ~ 150 Myr. For HR 8799c, the 3-cr ranges are T off = 
925 - 1025 K, log 10 c/ = 3.8 - 4.3, M = 6 - 13 Mj, and 
t = 20 - 100 Myr. The best-fit solution has T eff = 1000 
K, \og w g = 4.2, M ~ 11 Mj, and t ~ 65 Myr. Finally, 
for HR 8799d, the 3-cr constraints are T off = 850 - 975 
K, log 19 5 = 3.5 - 4.2, M = 3 - 11 Mj, and t = 10 - 70 
Myr, with the best fit obtained for T c g = 900 K, log 10 g 
= 3.8, M ~ 6 Mj, and t ~ 20 Myr. The best-fit ages 
for the three objects lie between 20 and 150 Myr. The 
age of the star has been reported to be 30 - 160 Myr 
(Marois et al. 2008), which is consistent with the age 
ranges of the objects, suggesting that the objects might 
have formed in situ, and that they are coeval. The new 



cloudy models developed in this work yield slightly im- 
proved fits over those previously reported by Curric ct 
al. (2011) for the same data, and our constraints on T e g 
and log 10 g, though slightly lower, are generally similar. 

We note that the observations at 3.3 /im for HR 8799b 
and HR 8799c indicate fluxes that are higher, by over 2-cr, 
than those predicted by the best-fitting models. This dif- 
ference can in principle be explained by a lower methane 
abundance than the equilibrium abundances considered 
in our cloudy models. 

Motivated by our results on the HR 8799 companions, 
we have provided a new suite of cloudy atmospheric mod- 
els in anticipation of the detection of similar objects in 
the future. These models are characterized by AE-type 
clouds that were able to fit the observations of HR 8799b, 
c, and d. We reported model spectra, photometric mag- 
nitudes in the infrared, and physical properties of the 
corresponding thick-cloud models. This model grid spans 
T cff = 600 - 1700 K, log lo5 = 3.5 - 5.0, M = 2 - 40 
Mj, and t = 4 Myr - 1 Gyr. The models assume so- 
lar abundances and chemical equilibrium for the gas- 
phase species. The models predict that objects like the 
HR 8799 planets will track a distinctly different locus in 
color- magnitude space from those of the L/T dwarf se- 
quence. The objects would diverge in color-magnitude 
space from the L/T dwarf sequence at the L/T tran- 
sition, and extend to systematically fainter and redder 
regions. We hope that the present models will help in- 
terpret the new observations of directly-imaged substcl- 
lar objects anticipated in the near future, to characterize 
their masses and ages, and constrain their planetary na- 
ture (Spiegel et al. 2011). The grid of models reported 
in this work, and those used in Currie et al. (2011), can 
be obtained in electronic format at the follow ing URL: 
http T/l www .astro . pr incet on . edu/ -burrows /[ 
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TABLE 1 

Physical parameters of best-fit models of the HR 8799 a planets 



Object 


T cB (K) 


l°Sio(3) b 


Age (Myr) 


Mass (Mj) 


Seperation (AU) C 


Model Name d 


HR 8799b 


~ 850 


~ 4.3 


~ 150 


~ 12 


68 


AE60. 850.1s. g4.3 


HR 8799c 


~ 1000 


~ 4.2 


~ 65 


~ 11 


38 


AE60.1000.1s.g4.2 


HR 8799d 


~ 900 


- 3.8 


~ 20 


~ 6 


21 


AE60. 900.1s. g3. 8 



a Stellar parameters: [Fe/H] = -0.47, Age = 30 - 160 Myr, Spectral type = A5V, Mass = 1.5 ± 0.3 M Q (Marois et al. 2008) 

g is in cm/s 2 . 
c Projected orbital separation (Marois ct al. 2008) 

^ The nomenclature 'AE60.850.1s.g4.3' denotes AE-typc cloud model, with 60-/im modal particle size, T e ff — 850 K, lxsolar mctallicity, and 
logio 9 = 4.3. 

° In cases where multiple solutions are available at similar levels of fit, a representative best-fit is chosen. 



TABLE 2 

Ranges of physical parameters constrained by the observations' 1 



Object T cB (K) log 10 (9) b Age (Myr) Mass (Mj) X 2 c 



HR 8799b 750-850 3.5-4.3 10-150 2-12 7.2-21.1 

HR 8799c 925-1025 3.8-4.3 20-100 6-13 14.2-21.0 

HR 8799d 850-975 3.5-4.2 10-70 3-11 14.9-23.5 



a The ranges for all the quantities reported here correspond to fits at 3-cr confidence, as described in section l4.ll and as shown in Fig. [6] Fig. ED 

and Fig. [TT] 

a g is in cm s -2 . 

c x 2 is denned in section l4.ll 



Ireland, M. J., Kraus, A., Martinache, F., Law, N., Hillenbrand, 

L. A. 2010, ApJ, in press, arXiv:1011.2201 
Janson, M., Bergfors, C, Goto, M., Brandner, W., & Lafrenire, 

D. 2010, ApJ, 710, L35 
Kalas, P., et al. 2008, Science, 322, 1345 

Lafreniere, D., Jayawardhana, R., & van Kerkwijk, M. H. 2008, 
ApJ, 689, L153 

Lafrenire, D., Marois, O, Doyon, R., & Barman, T. 2009, ApJL, 
694, L148 

Lafreniere, D., Jayawardhana, R., & van Kerkwijk, M. H. 2010, 

ApJ, 719, 497 
Lagrange, A.-M., et al. 2009, A&A, 493, L21 
Lcggett, S., et al. 2010, ApJ, 710, 1627 

Liu, M. O, et al. 2010, Proc. of the SPIE, Vol. 7736, 77361K 
Macintosh, B. et al. 2006, Proc. of the SPIE, Vol. 6272, 62720L 
Marley, M., Seager, S., Saumon, D., Lodders, K., Ackcrman, A. 
S., Freedman, R. S., & Fan, X. 2002, ApJ, 568, 335 



Marois, O, et al. 2008, Science, 322, 1348 

Marois, O, Zuckcrman, B., Konopacky, Q. M., Macintosh, B., & 

Barman, T. 2010, Nature, 468, 1080 
Metchev, S., & Hillenbrand, L. A. 2006, ApJ, 651, 1166 
Mohanty, S., Jayawardhana, R... Hulamo, N., & Mamajek, E. 

2007, ApJ, 657, 1064 
Saumon, D., Chabrier, G., & Van Horn, H.M. 1995, ApJS, 99, 713 
Saumon, D., Marley, M. S., Lodders, K., & Freedman, R. S. 2003, 

IAUS, 211, 345 
Saumon, D., & Marley, M. 2008, ApJ, 689, 1327 
Sharp, CM. & Burrows, A. 2007, ApJS, 168, 140 
Spiegel, D. S., Burrows, A., & Milsom, J. A. 2011, ApJ, 727, 57 
Sudarsky, D., Burrows, A., & Pinto, P. 2000, ApJ, 538, 885 
Tamura, M. 2009, AIP Conf. Series, 1158, 11 
Tsuji, T. 2002, ApJ, 575, 264 
Tsuji, T. 2005, ApJ, 621, 1033 



Model atmospheres for massive gas giants with thick clouds 11 



TABLE 3 

Physical and photometric characteristics of objects with thick "AE type" CLOUDS a ' b 



T eff (K) 


login 3 (cm 2 /s) 


Mass (Mj) 


Age 


(Gyr) 


z 


i 


Y 


J 


11 


Ks 


L' 


M 


600 


3 


5 


2.17 





,023 


21, 


74 


19 


.62 


IT 


83 


16. 


74 


1 5 


39 


13 


. 73 


12 


.25 


650 


3 


5 


2.23 





.017 


21 


11 


19 


. 1 1 


IT 


35 


16. 


14 


11 


85 


13 


.26 


12 


.02 


700 


3 


5 


2.30 





.014 


20 


5T 


18 


.58 


16 


84 


15. 


56 


1 1 


34 


12 


. 86 


1 1 


.85 


750 


3. 


5 


2.38 





Oil 


20 


15 


18 


.23 


16 


52 


15. 


1 1 


13 


94 


12 


.46 


1 1 


.Tl 


800 


3 


5 


2.45 





.009 


19 


TO 


IT 


.85 


16 


18 


1 4 . 


Tl 


1 3 


5T 


12 


. 1 1 


1 1 


.62 


850 


3 


5 


2.51 





.007 


19 


30 


IT 


.51 


15. 


88 


14. 


37 


1 3 


23 


11 


.80 


1 1 


r r 7 
Xi i 


900 


3. 


5 


2.58 





.006 


18 


9-°, 


IT 


18 


15 


60 


14. 


0T 


12 


91 


1 1 


. 54 


1 1 


.56 


950 


3 


5 


2.65 





,005 


18 


58 


16 


.89 


15. 


3 1 


13. 


80 


12 


60 


11 


.31 


11 


.56 


1000 


3 


5 


2.72 





.004 


18 


22 


16 


.57 


15. 


0T 


13. 


55 


12 


32 


11 


. 10 


1 1 


.57 


600 


1 





5.85 





,125 


22 


02 


19 


.87 


18. 


1 1 


IT. 


1 1 


16 


02 


14 


.10 


12 


.60 


650 


4 





5.98 





,098 


21 


18 


19 


. 14 


IT 


18 


16. 


16 


15 


45 


13 


. 1 1 


12 


.39 


700 


4 





6.11 





,078 


20 


85 


18 


.87 


IT. 


18 


15. 


98 


1 1 


92 


13 


.30 


12 


.18 


750 


1 





6.24 





,064 


19 


32 


IT 


.54 


16 


OT 


15. 


00 


11 


16 


12 


. 86 


12 


.01 


800 


4 





6.38 





,052 


19 


89 


18 


.04 


16. 


42 


15. 


10 


11 


09 


12 


. 63 


1 1 


.89 


850 


4 





6.52 





,043 


19 


45 


IT 


.66 


16. 


09 


14. 


T3 


1 3 


T4 


12 


. 33 


1 1 


.81 


900 


4 





6.67 





,035 


19 


06 


IT 


.32 


15 


T9 


1 4 . 


11 


13 


42 


12 


.06 


1 1 


.75 


950 


4 





6.83 





,029 


18 


T5 


IT 


.06 


15. 


56 


14. 


15 


1 3) 


12 


11 


.81 


1 1 


.73 


1000 


4 





6.98 





,025 


18 


13 


16 


.78 


15 


32 


13. 


90 


12 


84 


11 


.60 


11 


.T2 


1050 


4 





7.12 





,021 


18 


11 


16 


.51 


15. 


08 


13. 


6T 


12 


58 


1 1 


. 4 1 


1 1 


.73 


1100 


4 





7.25 





,019 


IT 


80 


16 


.23 


11 


85 


13. 


16 


12 


35 


11 


. 23 


1 1 


.75 


600 


4 


5 


14.46 


1 


,005 


22 


19 


19 


.98 


18 


28 


IT. 


12 


16 


86 


14 


.46 


12 


.96 


650 


4 


5 


14.69 





,787 


21 


63 


19 


.53 


IT. 


8T 


16. 


89 


16 


22 


14 


.08 


12 


. i 4 


700 


4 


5 


14.93 





,619 


20 


91 


18. 


.90 


IT 


32 


16. 


33 


15 


65 


13 


. i 3 


12 


. 54 


750 


4 


5 


15.15 





,510 


20 


11 


18 


.50 


16 


96 


15. 


89 


15 


16 


13 


.42 


12 


.3 / 


800 


4 


5 


15.38 





,434 


20 


IT 


18 


.31 


16. 


75 


15. 


52 


14. 


69 


13 


.09 


12 


.20 


850 


4 


5 


15.60 





,373 


19 


60 


IT 


.80 


16. 


30 


15. 


12 


11 


34 


12 


. 84 


12 


.10 


900 


4 


5 


15.83 





,324 


19 


23 


IT 


19 


16 


02 


14. 


T9 


11 


00 


12 


.58 


12 


.01 


950 


4 


5 


16.07 





,281 


18 


90 


IT 


.20 


15. 


77 


14. 


51 


13 


68 


12 


.34 


11 


.94 


1000 


4 


5 


16.33 





,245 


18 


59 


16 


.94 


15. 


53 


14. 


25 


1 3) 


40 


12 


. 12 


1 1 


.90 


1050 


4 


5 


16.57 





,216 


18 


30 


16 


.69 


15. 


31 


1 4 . 


02 


13 


13 


1 1 


.91 


1 1 


.89 


1100 


4 


5 


16.82 





,192 


18 


01 


16 


IT 


15. 


12 


13. 


82 


12 


8T 


11 


. i 3 


1 1 


.89 


1200 


4 


5 


17.37 





,149 


IT 


50 


16 


.00 


11 


Tl 


13. 


11 


12 


41 


11 


. 1 1 


1 1 


.92 


1300 


4 


5 


17.93 





117 


16 


9T 


15 


.54 


11. 


32 


13. 


1 1 


12 


10 


1 1 


. 10 


1 1 


.91 


1400 


4 


5 


18.42 





,096 


16 


45 


15 


.09 


13 


94 


12. 


T9 


11 


82 


10 


.85 


1 1 


.80 


1500 


4 


5 


18.95 





,079 


15 


89 


14 


.61 


13 


53 


12. 


•17 


11 


55 


10 


. 63 


11 


.65 


1600 


4 


5 


19.51 





,066 


15 


40 


14 


.21 


13 


IT 


12. 


18 


11 


32 


10 


.45 


1 1 


.46 


1700 


4 


5 


20.04 





,057 


11 


91 


13 


.83 


12 


83 


11. 


92 


11 


12 


10 


.30 


1 1 


.23 


600 


5 





32.78 


1 


,569 


22 


16 


20 


.11 


18. 


45 


IT. 


TT 


IT 


88 


14 


. 83 


13 


.36 


650 


5 





33.12 


3 


,627 


21 


Tl 


19 


.50 


IT 


92 


IT. 


18 


IT 


20 


14 


.4 1 


13 


.15 


700 


5 





33.45 


2 


938 


21 


31 


19 


.23 


IT. 


66 


16. 


T3 


16 


18 


14 


.09 


12 


.91 


750 


5 





33.80 


2 


,417 


21 


12 


19 


.04 


IT. 


50 


16. 


55 


16 


24 


13 


.94 


12 


.81 


800 


5. 





34.16 


2 


,021 


20 


32 


18. 


.37 


16 


90 


15. 


89 


15 


18 


13 


.51 


12 


.57 


850 


5. 





34.51 


1 


711 


19 


8T 


IT 


.99 


16. 


56 


15. 


53 


15 


0T 


13 


.30 


12 


.13 


900 


5. 





34.88 


1 


,463 


19 


48 


IT 


.68 


16 


28 


15. 


21 


11 


68 


13 


.06 


12 


.31 


950 


5. 





35.28 


1 


256 


19 


10 


IT 


.35 


15. 


99 


14. 


90 


11 


35 


12 


.85 


12 


.22 


1000 


5. 





35.70 


1 


078 


18 


80 


IT 


.12 


15. 


77 


14. 


64 


11 


03 


12 


.63 


12 


.14 


1050 


5. 





36.18 





,920 


18 


51 


16 


.87 


15. 


55 


14. 


40 


13 


Tl 


12 


.13 


12 


.09 


1100 


5. 





36.73 





768 


18 


23 


16 


.64 


15. 


31 


14. 


18 


13 


48 


12 


.25 


12 


.06 


1200 


5. 





37.54 





583 


IT 


T3 


16 


.21 


11. 


9T 


13. 


T9 


13 


01 


11 


.91 


12 


.05 


1300 


5. 





38.04 





,503 


IT 


2T 


15 


.83 


11. 


63 


13. 


IT 


12 


61 


11 


.62 


12 


.09 


1400 


5. 





38.62 





442 


16 


T9 


15 


.12 


11 


28 


13. 


IT 


12 


29 


11 


.31 


12 


,0T 


1500 


5. 





39.70 





360 


16 


31 


15 


.03 


13 


93 


12. 


89 


12 


03 


11 


.10 


11 


.97 


1600 


5. 





40.84 





,294 


15 


83 


14 


.62 


13 


58 


12. 


61 


11 


T8 


10 


.90 


11 


.84 


1700 


5. 





41.62 





260 


15 


10 


14 


.27 


13. 


25 


12. 


35 


11 


5T 


10 


.Tl 


11 


.66 



a All the models are assumed to have solar metallicity. 

b Absolute photometric magnitudes arc reported in the following filters: z' (SDSS), Y(MKO), J (MKO), H (MKO), Ks (ISAAC), L' (MKO), and 
M (MKO). 
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Fig. 1. — Model spectra of giant planet atmospheres without clouds. Th e so lid lines show model spectra, and photometric observations 
of HR 8799b are shown in black. The observations are described in section [47T] Each panel shows a range of models exploring a particular 
parameter: T e ff , log 10 (g), [Fc/H], and the eddy diffusion coefficient K zz . The models in the top-left, top-right, and bottom-left panels assume 
equilibrium chemist ry w ith the specified abundances. The models in the bottom-right panel have non-equilibrium CO-CH4 chemistry, as 
described in section f3.f I and in Hubeny & Burrows (2007). The model names in the legend specify the parameters values. For example, in 
the top-left panel, T800.1s.g4.5 denotes a cloud-free ('T') model with T c g = 800 K, log 10 (g) = 4.5, and f xsolar metallicity. For the models 
in the bottom-right panel, the abundances are solar, and the legend specifies the diffusion coefficient; for example, 'd6' implies K zz = fO 6 
cm 2 s _1 . Cloud-free models predict ~2-5 times higher fluxes in the z, J, and H bands than are observed for HR 8799b. All three objects 
considered in this study, HR 8799b, c, and d, require thick clouds to explain the observations. 
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1234567 1234567 
X(/xm) X(ju,m) 



Fig. 2. — Model spectra of giant planet atmospheres with clouds. Four types of cloudy models, as described in section [2] are shown. 
The top-left panel shows models with E-type clouds which truncate at high altitudes. The top-right panel shows models with A-type 
clouds which extend up to the top of the atmosphere, and have the thickest extent among the models considered in this work. The bottom 
panels show AEE and AE cloud types, which are intermediate between the A and E types, as described in section [2] The nomenclature 
E60.800.1s.g4.0 stands for E-type clouds with 60-/^m modal particle size, T e ff = 800 K, lxsolar metallicity, and log 10 g = 4.0. Obs ervat ions 
of HR 8799b arc shown in black. AE-type clouds provide the best match to the observations (see Fig. [4] and discussion in section [4.2l l. 
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500 1000 1500 2000 2500 3000 3500 

T(K) 

Fig. 3. — Temperature-Pressure (T-P) profiles of model atmospheres with clouds. The model nomenclature shown in the legend is 
explained in Fig.[2]and in section[2] All the models have T c ff = 1000 K, log 10 g = 4.0, and solar metallicity, but different cloud types. The 
A-type (E-type) forsterite clouds which have the thickest (thinnest) cloud cover have the hottest (coolest) T-P profiles for the same T e g . 
Models with AE-type forsterite clouds with 60-/^m modal particle size and those with A-type iron clouds with 1-pm modal particle size, 
with 1% supersaturation, have similar T-P profiles. 
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Fig. 4. — Model spectra of HR 8799b. The obs ervat ions, shown in black, arc from the literature (described in section [4. 111 . The model 
nomenclature is described in Fig. [2] and in section [4. 21 The black lines at the top of each panel indicate the photometric bandpasses, and 
the colored solid lines show the model spectra. The top-left panel shows the best-fit models. The observations are best explained by cloudy 
models with intermediate cloud cover (AE type), with modal particle size of 60 /an, T c g = 850 K, metallicity of lx solar, and log 10 (g) = 4.3, 
shown in the red curve. The blue filled circles show the best-fit model integrated over the photometric bandpasses. The brown curve in 
the top-left panel shows a model with Fe clouds with 1% supersaturation, and 1-^tm particle size, with T e ff = 800 K, metallicity of lx 
solar, and log 10 (g) = 4.0, which also provides a good fit to the observations, showing that different cloud types and compositions can in 
principle explain the photometric points. The remaining panels show departures from the best-fit model (red model in top-left panel) in 
temperature (top-right), gravity (bottom-left), and cloud type (bottom-right). The values of the varied parameter for the different models 
are shown in the legend; all unspecified parameters are the same as for the best-fit model. Higher T e g yields systematically higher fluxes in 
all the channels, except at 3.3 fim. E-type cloud models which explain the spectra of L and T dwarfs fairly well are inadequate to explain 
the current observations. For a given T e g, E-type cloud models predict 2-5 times higher fluxes at shorter wavelengths (in z, J, and H 
bands), than those required by the data. On the other hand, purely A models yield excess flux at short wavelengths and low flux at longer 
wavelengths. The intermediate AE-type models explain the data best. Changing the metallicity does not improve the fits significantly. See 
text (section 14.21 1 for discussion. 
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600 800 1000 1200 1400 10 20 30 40 

T eff (K) Mass(M J ) 

Fig. 5. — Model space explored for HR 8799b. The symbols show models with AE clouds computed over a grid in T e ff and log 1 g(g), as 
shown in the left panel. The right panel shows the corresponding regions in the Mass- Age plane, derived from T c ff and log 1C )(g) (Burrows 
et al. 2001). The models are color coded by the degree of fit to the data, quantified by the confidence level associated with the \ 2 ( see 
section [4.11 1. The x 2 i s derived from fits to photometric data in all the nine wavebands available. The purple open-diamonds show fits 
to within the 3cr confidence level, and the red, blue, green, black, and gray fillcd-circles indicate fits within 4<r, 5<r, 6<r, 7a, and > 7cr, 
confidence levels, respectively. The colored points enclosed in black boxes indicate models that yield an age for the object that is greater 
than the reported maximum age of the star (160 Myr), and, hence, are ruled out as implausible solutions despite their level of fit to the 
data. 
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Fig. 6. — Optimal constraints on the atmospheric parameters of HR 8799b. The description of this figure is same as that for Fig. [5] with 
the exception that the 3.3-/mi observation is excluded from the x 2 i implying 7V (, S = 8 observations and Ndof = 7 degrees of freedom (sec 
section |32] for a discussion). 
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Fig. 7. — Model spectra of HR 8799c. The observat ions are shown in black. The top-left panel shows the best fitting models. The 
model nomenclature is described in Fig.[2]and in section ^, 21 AE-typc clouds with 60-fim modal particle size (sec section[2]| provide a good 
match to the observations, similar to the case of HR 8799b shown in Fig. [4] The blue filled circles show the AE model integrated over the 
photometric bandpasses. Iron clouds, with 1% supcrsaturation and l-/im modal particle size, also explain the observations fairly well. The 




best fitting AE models require T c g ~ 1000 K and 



~ 4.2; current data are not sensitive to changes in metallicity. The remaining 



three panels show the dependence of model spectra on the key atmospheric parameters, T e g , log 1 g(g), and cloud type. T c ff > 1 000 K leads 
to substantially higher 3.8-fiia flux than observed, whereas the J-band observation rules out models with lower T ef j. See section f4. 31 in text 
for details. 
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Fig. 8. — Model space explored for HR 8799c. See Fig.[5]for description of the symbols and colors. 
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Fig. 9. — Optimal constraints on the atmospheric parameters of HR 8799c. The symbols and colors arc described in Fig. [5] The parameter 
space explored is same as that in Fig. [H] Howev er, in th e con straints shown in this figure we exclude the contribution of the 3.3-^jm data 
point to the x 2 > for reasons discussed in sections 14.21 and 14.31 
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Fig. 10. — Models fits to observations of HR 8799d. The observations are shown in black. Three best fitting models with different cloud 
types are shown in the top-left panel. The model nomenclature is described in Fig. [2]and in section|4]2] The observations can be explained 
by a wider range of models than those for HR 8799 b and c. Nevertheless, thick clouds are required to explain the observations. As shown 
in the top-left panel, models with A-type clouds as well as AE-type clouds provide almost equally good fits to the data. The red curve in 
the top-left panel shows a best-fit model with AE-type cover, with T e ff ~ 900 K and log 10 (g) ~ 3.8; the blue filled circles show the flux 
integrated in the photometric bandpasscs. Iron clouds, with 1-fim modal particle size and 1% supcrsaturation, also provide a good match 
to the observations, as observed for HR 8799 b and c. The remaining panels show the dependence of the spectra on different atmospheric 
parameters. 
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Fig. 11. — Optimal constraints on the atmospheric parameters of HR 8799d. The approach is described in Fig. [5] and Fig- [6j For 
HR 8799d, inclusion of the 3.3 /im photometric point in the \ 2 does not significantly alter the x 2 , since only an upper limit was observed 
in this bandpass which is higher than all the models explored. 
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Fig. 12. — Color-magnitude diagrams comparing the loci for standard brown dwarf atmosphere models with cloud decks (Case "E") and 
loci for planet-mass objects with thick clouds (Case "AE"). Overplotted are data for HR 8799 (Currie et al. 2011), other planet-mass 
objects, and field L/T dwarfs from Leggett et al. (2010), shown in grey filled circles. In all diagrams, the thick cloud models predict a far 
different distribution of colors and magnitudes for cool, substellar objects. We predict that many 1-10 Mj planets that will be discovered 
with GPI and SPHERE will occupy regions of color-magnitude space coincident with the thick cloud model locus and define a "planet 
locus" distinguishable from the L/T dwarf sequence. 



